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Abstract — A  new  frequency-domain  channel  estimation  and 
equalization  (FDE)  scheme  is  proposed  for  single  carrier  (SC) 
underwater  acoustic  communications.  The  proposed  SC-FDE 
employs  a  small  training  signal  block  for  initial  channel  es¬ 
timation  in  the  frequency  domain  and  converts  the  estimated 
transfer  function  to  a  desired  DFT  (discrete  Fourier  transform) 
size  for  channel  equalization  of  the  data  blocks.  The  frequency 
domain  equalizer  is  designed  using  the  linear  minimum  mean 
square  error  criterion.  A  new  phase  coherent  detection  scheme 
is  also  proposed  and  deployed  to  combat  the  phase  drift  due 
to  the  instantaneous  Doppler  in  the  underwater  channels.  The 
channel  transfer  functions  and  group-averaged  phase  drift  are 
re-estimated  adaptively  in  a  decision-directed  manner  for  each 
data  block  in  a  packet,  which  contains  M  blocks  of  QPSK  data. 
The  proposed  SC-FDE  method  is  applied  to  single  input  multiple 
output  (SIMO)  systems  using  the  experimental  data  measured 
off  the  coast  of  Panama  City,  Florida,  USA,  June  2007.  The 
uncoded  bit  error  rate  of  the  SIMO  systems  varies  between  1.3% 
to  6.8  x  10  5  when  4  ~  8  receive  hydrophones  are  utilized,  and 
the  source-receiver  range  is  5.06  km. 

I.  Introduction 

Shallow  water  horizontal  communication  channels  are  often 
hostile  for  high  data  rate  underwater  communications  due  to 
two  major  obstacles.  One  is  the  excessive  multipath  delay 
spread  in  a  medium  range  shallow  water  channel  which  is 
usually  on  the  order  of  10-30  ms  and  causes  the  intersymbol 
interference  (ISI)  to  extend  over  20-300  symbols  at  a  data 
rate  of  2-10  kilosymbols  per  second.  Another  obstacle  is  the 
time-varying  Doppler  shift  due  to  relative  motion  between 
the  source  (transducers)  and  receiver  (hydrophones),  dynamic 
motion  of  water  mass,  and  varying  sound  speed,  etc.  A  few 
Hz  of  Doppler  in  underwater  channels  can  be  very  significant 
due  to  the  low  velocity  of  acoustic  waves  (~1500  m/s).  The 
ratio  of  Doppler  to  carrier  frequency  in  underwater  channels  is 
on  the  order  of  1 0  to  10-4;  while  the  ratio  in  RF  wireless 
channels  is  on  the  order  of  10-'  to  HD9.  The  significant 
Doppler  shift  cause  not  only  rapid  fluctuation  in  the  fading 
channel  response  but  also  compression  or  dilation  of  signal 
waveforms.  These  two  obstacles  make  the  coherent  receiver 
of  underwater  communication  systems  much  more  complex 
[l]-[4]  than  RF  systems. 

It  has  been  successfully  demonstrated  in  [1]  that  coherent 
detection  of  underwater  acoustic  communications  can  be  real¬ 
ized  by  joint  decision  feedback  equalization  (DFE)  and  phase 
synchronization  that  employs  a  phase-locked  loop  (PLL)  or 
delay-locked  loop  (DLL).  However,  the  DFE  and  PLL/DLL 


have  to  interact  in  a  nonlinear  fashion  and  in  a  symbol- 
by-symbol  basis,  therefore  it  requires  careful  selection  of 
the  number  of  equalizer  taps  and  tuning  of  the  equalizer 
and  PLL/DLL  coefficients.  Stable  and  robust  operation  of 
the  time-domain  equalizer  is  sometimes  difficult  to  obtain  in 
different  channel  conditions.  Recent  improvement  on  robust 
time-domain  DFE  has  been  reported  in  [4]  using  a  fixed  set  of 
parameters  at  a  cost  of  slightly  degraded  bit  error  performance. 

Further  increasing  of  data  rates  imposes  a  greater  challenge 
on  time-domain  coherence  receivers.  Therefore,  orthogonal 
frequency  division  multiplex  (OFDM)  techniques  have  been 
applied  to  underwater  acoustic  channels  [5],  [6]  with  some 
promising  results.  Adaptive  channel  and  phase  estimation  was 
used  in  [5]  to  update  the  channel  Doppler  and  carrier  frequency 
offset  (CFO)  from  block  to  block  and  secure  successful  data 
detection.  In  [6],  carrier  frequency  offset  and  channel  estima¬ 
tion  were  performed  within  each  OFDM  block  employing  25% 
of  the  OFDM  bandwidth  for  pilot  tones  and  75%  for  payload. 
Therefore  it  is  suitable  for  fast  varying  channels  provided 
that  the  frequency  selective  channel  length  is  less  than  the 
number  of  pilot  tones.  Both  algorithms  employ  quadrature 
phase  shift  keying  (QPSK)  and  assume  constant  Doppler 
shift  and/or  carrier  frequency  offset  (CFO)  over  one  OFDM 
block.  It  has  been  shown  by  field  data  measurements  that  the 
OFDM  schemes  are  successful  for  block  sizes  up  to  1024  data 
symbols. 

Alternative  to  the  OFDM  approach,  a  single  carrier  sys¬ 
tem  with  frequency-domain  equalization  (SC-FDE)  has  been 
shown  [7],  [8]  to  have  comparable  performance  and  com¬ 
plexity  for  RF  broadband  wireless  communications.  Besides, 
SC-FDE  has  the  advantage  of  lower  peak-to-average  power 
ratio  and  less  sensitivity  to  carrier  frequency  offsets.  Compared 
to  time-domain  equalization,  SC-FDE  has  better  convergence 
properties  [9]  and  less  computational  complexity  in  severe 
inter-symbol  interferences  (ISI)  channels.  However,  the  SC- 
FDE  schemes  developed  for  broadband  RF  uplink  commu¬ 
nication  systems  may  not  be  applied  directly  to  underwater 
channels. 

In  this  paper,  we  propose  a  new  frequency-domain  channel 
estimation  and  equalization  scheme  for  single  carrier  underwa¬ 
ter  acoustic  communications.  The  proposed  SC-FDE  employs 
a  small  training  signal  block  for  initial  channel  estimation 
in  the  frequency  domain  and  converts  the  estimated  transfer 
function  to  a  desired  DFT  (discrete  Fourier  transform)  size  for 
channel  equalization  of  the  data  blocks.  The  frequency  domain 
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equalizer  is  designed  using  the  linear  Minimum  Mean  Square 
Error  (MMSE)  criterion.  A  new  phase  coherent  detection 
scheme  is  then  developed  to  combat  the  phase  drift  due  to 
the  instantaneous  Doppler  in  the  underwater  channels.  The 
channel  transfer  functions  and  group-averaged  phase  drift  are 
re-estimated  adaptively  in  a  decision-directed  manner  for  each 
data  block  of  a  packet,  which  contains  M  blocks  of  QPSK 
data.  The  proposed  SC-FDE  method  is  applied  to  SIMO 
systems  using  the  experimental  data  measured  off  the  coast 
at  Panama  City,  Florida,  USA  in  June  2007.  The  uncoded 
bit  error  rate  of  the  SIMO  systems  varies  between  1.3%  to 
6.8  x  10-5  when  4  ~  8  receive  hydrophones  are  utilized,  and 
the  source-receive  range  is  5.06  km. 

The  advantages  of  the  proposed  SC-FDE  scheme  are  two 
folds.  First,  the  size  of  the  training  signal  block  can  be  different 
from  the  DFT  size.  As  long  as  it  covers  the  fading  channel 
length,  it  can  be  selected  small  to  reduce  the  overhead.  Second, 
the  new  phase  correction  scheme  operates  on  averaged-group 
phase  drift  rather  than  individual  symbol’s  phase  thus  provid¬ 
ing  robustness  against  noise  and  interference. 

This  paper  also  presents  a  generic  channel  model  for  SISO 
and  SIMO  underwater  acoustic  communication  systems.  A 
rigorous  derivation  of  the  proposed  SC-FDE  is  then  developed 
based  on  the  model  and  the  mathematical  formulas  provide 
theoretical  support  and  insightful  understanding  on  how  the 
method  achieves  the  high  performance. 

II.  The  System  Model  and  Preliminaries 

Consider  a  underwater  acoustic  communications  system 
using  a  single  transducer  source  and  M -hydrophone  receiver. 
The  baseband  equivalent  signal  received  at  the  m-th  hy¬ 
drophone  can  be  expressed  in  discrete-time  domain  as 

L 

ym{k)  =  +vm(fc)  (1) 

1=1 

where 

k  time  index; 

ym{k)  received  signal  sample; 

x(k)  transmitted  data  symbol  or  pilot  symbol; 

hm(l,k)  impulse  response  of  the  frequency-selective, 

time-varying  fading  channel  with  length  L; 

L  fading  channel  length  in  terms  of  T; 

T  data  symbol  interval; 

fm  k  time-varying  carrier  frequency  offset  (CFO)  or 

instantaneous  Doppler  shift  caused  by  relative 
motion  between  the  transducer  and  hydrophone, 
dynamic  motion  of  water  mass,  varying  sound 
speed,  and  changing  water  columns,  etc.; 

6m,o  phase  error  after  coarse  synchronization; 

Vm{k )  white  Gaussian  noise  with  average  power  er2. 

If  the  Doppler  shift  fm  k  is  significant,  then  it  causes  the 
received  signal  ym(k)  to  be  time-scaled  (compressed  or  di¬ 
lated)  signal  [2],  [10].  In  this  case,  re-scaling  and  re-sampling 
are  required  before  equalization  taking  place  to  cancel  1ST 


If  the  relative  motion  between  the  source  and  receiver  is 
insignificant,  then  the  average  Doppler  shift  E{fm  fc}  is  close 
to  zero,  but  the  instantaneous  Doppler  fmk  is  a  time-varying 
random  variable  with  a  zero  mean.  In  this  paper,  we  focus 
on  the  case  where  the  source  and  receiver  are  stationary  with 
E{fm  fc}  =  0.  The  case  with  E{fm  0}  ^  0  will  be  considered 
in  our  future  work. 

It  is  noted  that  the  fading  coefficients  hm(l,k )  include  the 
effects  of  the  transmit  pulse-shape  filter,  physical  multipath 
fading  channel  response,  and  the  receive  matched  filter.  It 
is  found  that  the  fading  channel  coefficients  usually  change 
much  slower  than  the  instantaneous  phase  27t  fm  k  kT  in  many 
practical  underwater  acoustic  channels  [1],  [3],  That  is  the 
main  reason  to  explicitly  separate  the  phase  27r  fm  k  kT  from 
the  fading  coefficients  hm(l,k)  in  (1). 

To  facilitate  frequency-domain  channel  estimation  and 
equalization  for  the  system  described  in  (1),  the  transmitted 
data  sequence  {x(fc)}  are  partitioned  into  blocks  of  size  N. 
Each  block  is  then  added  with  Nzp  cyclic  prefix  or  padded  with 
Nzp  zeros  [9],  For  the  zero-padding  scheme,  a  transmission 
block  is 

xbk  =  [  x(l)  ar(2)  •••  x(N)  0  •••  0  ]‘  (2) 

The  corresponding  received  signal  ybk  of  the  block  is  denoted 
by 

ybm  =  [y(1)  2/(2)  •••  y(N)  •••  viN+N^y  (3) 

where  the  superscript  [■]*  is  the  transpose.  The  length  of  zero 
padding  Nzp  is  chosen  to  be  at  least  L  —  1  so  that  the  inter¬ 
block  interference  is  avoided  in  the  received  signal. 

Adopting  the  overlap-add  method  [11]  for  iV-point  DFT,  we 
define  two  signal  vectors 

x  =  [  x(l)  x(2)  •••  x(N)  Y  (4) 

ym  =  [l/m(l)  ymi^Ezp)  t/m  (-^2p-T  1)  '''  ] 

+  [ym(N+l)  ■■■  ym(N+Nv)  0  •••<)]*  (5) 

and  their  corresponding  frequency-domain  signals  Ym  = 
FNym  and  X  =  Fv x,  where  Fv  is  the  normalized  DFT 
matrix  of  size  N  x  N,  i.e.,  its  (to,  n)-th  element  is  given  by 
exp  .  Note  that  F^FW  =  I*. 

The  time-domain  signals  x  and  ym  are  related  as 

y m  =  TmDmx  -T  vm  (6) 

with 

-MM)  o  •••  o  MM  •••  M2,i)- 
M 2,2)  MM)  o 

Tm=hm(L,L)  MM)  0  \  0 

0  h>m(L,L-\- 1)  ■  -  ... 

:  0 
-  o  •  •  •  0  hJL,N)  .  MMV) . 

(V) 
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Dm=diag|ei(27r/-.iT^m’o)  •••  ej{2nf™^NT'¥,m’o)}  (8) 

Vm  =  [l?m(l)  **■  Wri.  (Xzpj  nm(AAp- t”l)  ***  ^m(-^)] 

+  [um(W+l)  •••  vm(N+Nv)  0  •••  0]*.  (9) 

The  frequency-domain  representation  is 

Y  m  =  FNTmFjFNDmFjFlvx  +  FNvm 

=  Hm$roX  +  Vm  (10) 

where  Hm  =  Fv  T„,F^'  is  the  frequency-domain  channel 
matrix  (or  channel  transfer  function  matrix)  for  the  block,  and 
d>m  =  Fa,  DmF^  the  frequency-domain  Doppler  and  phase 
matrix  for  the  block. 

Although  <t>m  is  generally  a  non-diagonal  matrix,  its  diag¬ 
onal  elements,  n)}^=1  are  identical  and  equal  to 

®m{n,ri)=  — ^e-7(27r/"*.fcfc'w",’o),  n=  1,2,---  ,7V.  (11) 

fc= t 

This  property  of  the  single-carrier  frequency-domain  represen¬ 
tation  enables  effective  channel  estimation  and  equalization. 

III.  Frequency-Domain  Channel  Estimation 
In  the  channel  estimation  mode,  a  small  block  of  Nts 
training  symbols  are  padded  with  Nzp  zeros  (N^  <  7Vte) 
and  transmitted  over  the  acoustic  channel.  We  choose  the 
number  of  training  symbols  Nts  such  that  the  training  block 
time  duration  (7Vte  +  Nzp)T  is  smaller  than  a  third  of  the 
channel  coherence  time.  Therefore,  the  non-diagonal  ele¬ 
ments  of  d>m,  i.e.,  with  n  l,  are  negligi¬ 

ble  comparing  to  the  diagonal  elements  <3>TO(n,  n).  More¬ 
over,  the  fading  channel  coefficients  remain  constant  for 
the  entire  training  block,  thus  the  time -domain  channel  ma¬ 
trix  Tm  becomes  a  circulant  matrix  and  the  frequency- 
domain  channel  matrix  Hm  becomes  a  diagonal  matrix. 
That  is,  Hm=diag  (iTm(l),  Hm( 2),  •••,  fTm(7Vfe)}  with 

Hm(n )  =  l  hm{l,  1)  exp  •  Therefore,  the 

frequency-domain  representation  (10)  can  be  simplified  as 

Ym  (n)  =Hm  (n)  $m(n,n)X  (n)  +  Vm  (n) 

=\mHm(n)X(n)  +  Vm(n),  n=  1,2,  ,NP(12) 

where  Am  =  l  eY2’nfm-kT+em’0\  is  a  complex-valued 

unknown  parameter  with  amplitude  close  to  unit. 

The  channel  transfer  function  A mHm(n)  is  estimated  by  the 
least  squares  criterion  as 

Y  (n) 

A mHm(n)  =  ^YX,  n  =  1, 2,  •  •  *  ,  7Vfe.  (13) 
X[n) 

The  estimate  XmHm(n)  can  be  further  improved  by  a 
frequency-domain  filter  to  reduce  noise 

AmHm  =  AroFJVfa  [1  :  Ay  (F^  [1  :  Hm  (14) 

where  Fw  [1  :  Ay  denotes  the  first  Nzp  columns  of  FN  , 
Hm  =  diag  {  Hm(l)  Hm( 2)  •••  Hm{Nts)  },  and  Hm  = 

diag  { Hm(l)  Hm{ 2)  ^m(7Vfe)}. 


The  Nis -point  transfer  function  AmHm  can  be  directly 
transformed  into  an  Appoint  transfer  function  by 

AmHm  =  AmFw  [1  :  Ay  (F^  [1  :  Ay)*  Hm  (15) 

where  Fw  [1  :  Ay  denotes  the  first  columns 

of  Fv  (the  Appoint  DFT  matrix)  and  Hm  = 
diag{irm(l)  7Tm(2)  •••  Hm(N)  }.  The  estimated 

A’-point  transfer  function  is  then  used  for  frequency-domain 
channel  equalization  which  is  detailed  in  the  next  section. 


IV.  Frequency-Domain  Channel  Equalization 

In  the  data  transmission  mode,  we  choose  the  data  block 
time  duration  (N + N^,)T  slightly  smaller  than  the  channel  co¬ 
herence  time,  therefore  the  frequency-domain  channel  matrices 
of  all  M  hydrophones  {H m}^f=1  are  diagonal  matrices.  Based 
on  (10),  a  block  of  the  received  signals  at  the  M  hydrophones 
can  be  expressed  in  the  frequency  domain  as 


'  Yx  - 

■  Hxd?!  ‘ 

‘  V  -,  ■ 

y2 

= 

H2d>2 

x+ 

V2 

.  ym  _ 

1 

§ 

is 

K 

_ 1 

_VM. 

The  training  block  is  utilized  to  estimate  the  channel  transfer 
function  matrices  {AmHm}^£=1  and  these  are  used  to  equalize 
the  first  data  block.  Applying  the  MMSE  criterion,  we  obtain 
the  frequency-domain  equalized  block  data  as 


x=  £|Amf  Hj,Hm  +<J2Ijv  J  KrMXrn)  +  V 

\m=l  /  \ra=l  / 

/  M  \  —  1  /  M  \ 

=  I  A™ .  1 : 2  HlHm  +a2iy  A^H5,H,„$m  X  +  V 


,771  =  1 

M 


=  5Z  Arn<t>rn  X  +  V 


(17) 


where 


(A;nHrftnHm)  (18) 


is  a  diagonal  matrix  due  to  the  diagonal  properties  of 

{Hm}m=t  and  {Hm}"=1. 

Applying  inverse  DFT  to  the  equalized  data  vector  X,  we 
have  the  time-domain  data  vector  x  given  by 


M 

*=  F^X=^F^Am$mX  +  F^V 

m= 1 
M 

=  ^F^Am(FNDmF^)(FJvx)  +  v 

771=1 

M 

=  Y,  (F^AmFJV)Dmx  +  v.  (19) 

771=1 


Since  Am  is  a  diagonal  matrix,  all  the  diagonal  elements  of 
(F^AmFN)  are  identical  and  equal  to  7 m  =  -^trace(Am). 
When  the  data  block  time  duration  is  less  than  the  channel 
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coherence  time,  all  the  non-diagonal  elements  of  (F^AmFw) 
are  insignificant  comparing  to  jm.  Therefore,  the  A'-th  symbol 
of  x  can  be  expressed  by 


x(k)  = 


M 

J2  lmA2^kT+0m,o) 


jn—1 


x{k)  +  v(k) 


|  e3Zl3kx(k)  +  v(k) 


(20) 


where  (3k  =  Em=i  1me3^fm’kkT+9m’°\ 

From  (20),  we  can  conclude  that  the  complex-valued 
symbol-wise  scaling  factor  (3k  is  actually  a  diversity  combining 
factor  determined  by  the  M  channel  transfer  functions,  time- 
varying  Doppler  and  timing-error  phases.  In  other  words,  the 
equalized  data  symbol  x(k)  is  an  amplitude-scaled  and  phase- 
rotated  version  of  the  transmitted  data  symbol  x(k).  The 
rotating  phase  Z/3k  is  a  collection  of  all  the  contributions  from 
the  instantaneous  Doppler  fm  k  and  timing-error  phases  9m  0  of 
all  the  M  fading  channels.  For  each  individual  fading  channel, 
the  rotating  phase  Zp)mik  =  2nfm  kkT  +  9m0  +  Z-ym,  which 
represents  the  m-th  channel’s  Doppler-driven  shifting  phase, 
timing-error  phase,  and  the  channel  transfer  function  effect. 
This  is  certainly  a  clear  physical  interpretation  for  the  single 
carrier  frequency-domain  equalized  data. 

To  equalize  the  next  block  of  data,  the  channel  transfer 
function  matrices  need  to  be  re-estimated  using  the  detected 
data  and  received  signals  of  the  current  block,  then  employ 
these  estimated  transfer  function  matrices  to  equalize  the  next 
data  block.  This  scheme  is  commonly  referred  to  as  decision- 
directed  scheme.  Details  are  omitted  here  for  brevity. 

If  x(k)  is  phase  shift  keying  (PSK)  modulated  data,  then 
the  time-varying  rotating  phase  Z(3k  must  be  compensated  at 
the  receiver  after  FDE  and  before  demodulation  and  detection. 
This  is  discussed  in  detail  in  the  next  section. 


V.  Phase-Coherent  Detection 

In  this  section,  we  present  a  new  algorithm  for  estimating 
the  phase  Z(3k,  which  is  crucial  for  successful  data  detection 
of  PSK  modulated  symbols.  The  challenge  of  this  phase 
estimation  is  that  we  face  to  M  fading  channels,  in  which 
each  individual  channel  has  different  timing-error  phase  and 
time-varying  Doppler,  and  the  rotating  phase  Zf3k  represents  a 
nonlinearly  composed  effect  of  these  random  (or  time-varying) 
factors  of  all  the  M  fading  channels.  Therefore,  directly 
estimating  these  Doppler  and  timing-error  phases  will  be  very 
costly  if  any  possible. 

In  the  literature  of  underwater  acoustic  communications, 
phase  tracking  is  commonly  carried  out  by  utilizing  first- 
order  or  second-order  phase-locked  loop  or  delay-locked  loop 
[1],[3],  and  it  is  often  jointly  done  with  decision  feedback 
equalizers.  However,  this  approach  is  usually  sensitive  to 
channel  conditions  [4]  and  noise. 

What  we  know  from  the  nature  of  ocean  waters  is  that  the 
instantaneous  Doppler  fm  k  changes  gradually  from  time  to 
time,  i.e.,  it  does  not  change  arbitrary  in  a  short  period  of  time. 
Therefore,  the  rotating  phase  Z/3k  is  also  changing  gradually 


from  time  to  time.  We  treat  Z/3k  to  be  a  constant  for  a  small 
number  of  Ns  consecutive  received  symbols,  and  to  another 
constant  in  the  next  Ns  consecutive  received  symbols,  and  so 
on  so  forth. 

We  partition  the  equalized  A'-symbol  block  data  x  into  Ng 
groups,  each  with  Ns  data  symbols,  except  for  the  last  group 
might  have  less  than  Ns  symbols  if  N/Ng  is  not  an  integer. 

Let  ipp  denote  the  estimated  aver¬ 
age  rotating  phase  for  the  p-th  group  of 
{  Z/3(p_i)jvs+i,  Zf3(p_ i)jvs+2,  •••)  ZLf3(p-i)Ns+Ns }.  with 
p  =  1,  2,  •  •  •  ,Ng,  and  let  ip o  denote  the  initial  rotating  phase, 
A ipp  the  phase  difference  ipp—ipp-\.  Hence 

ipp  =  ipp-i  +  Aipp,  p=  1,2,  •••  ,Ng.  (21) 


For  MPSK  modulation  with  symbols  taken  from  an  M- ary 


constellation  AM  = 


leXP  l  W  J  >  TO  =  1>  2.  •  •  •  .  Mr 

we  define  a  phase  quantization  function  Q  [<p]  as  follows 


(m  —  1)27t 

MZ  ’ 


mfl'K—Z'K  ^  i  s'  m2n-^K 

Mm  ^  V  -  Mm 
m=  1,2,  •••  ,Mm. 


(22) 


We  are  now  in  a  position  to  present  our  algorithm  as  follows: 

Algorithm:  Group-wise  Rotating  Phase  Estimation 

Step  1.  Designate  the  first  Np  symbols  of  each 

transmitted  block  data  x  as  pilot  symbols  for  phase  reference 
to  determine  the  initial  rotating  phase  ip0  given  by 

1  Np 

=  W  )  -  Zx(k).  (23) 

jvp  fc= l 

Compensate  the  phase  of  the  first  group  data  by  e~3^ot 
yielding 

xk  (k)  =  x(k)e~3^° ,  k  =  1, 2,  •  •  •  ,  Ns.  (24) 


Calculate  the  individual  phase  deviation  from  its  nominal 
phase  of  each  symbol  in  the  first  group 

<Pi,k  =  Zxk  (k)  -  Q  [Zxk  (fc)] ,  k  =  1,2,-  •  •  ,NS.  (25) 

Calculate  the  average  phase  deviation  and  compute  the 
rotating  phase  for  the  first  group  as  follows 


II 

<1 

1  Na 

(26) 

II 

s  k=  1 

ipo  +  Aipi. 

(27) 

Set  p  =  2  for  next  step. 


Step  2.  Compensate  the  phase  of  the  p-th  group  data  by 
e-ii>p-it  yielding 

xp(k)  =x((p-l)Ns+k)e~j^\  k=  1,2,  -  -  ,NS.  (28) 

Calculate  the  individual  phase  deviation  from  its  nominal 
phase  of  each  symbol  in  the  p-th  group 

<pPik  =  Zxp(k)  -  Q[Zxp(k)] ,  k=  1,2, •••  ,NS.  (29) 
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Calculate  the  average  phase  deviation  and  estimate  the 
rotating  phase  for  the  p- th  group  as  below 


II 

<3 

1  N‘ 

N~s  ^  ^ 

(30) 

V’p  = 

s  fc= i 

4>p- 1  +  A  ij}p. 

(31) 

Step  3.  Update  p  =  p  +  1,  repeat  Step  2  until  p  =  Ng.  □ 

After  estimating  the  Ng  group  phases,  we  can  compensate 
the  phase  rotation  of  the  equalized  data  x  in  group  basis: 

xp(k)  =  x((p-l)Ns+k)e~^,  pZit2,"'.’,Ng  •  O2) 

Finally,  the  binary  information  data  of  the  block  can  be 
obtained  via  standard  MPSK  demodulation  procedure  on  the 
phase-compensated  signal  xp(k)  of  the  block.  Next  block  data 
can  be  processed  in  a  similar  manner,  details  are  omitted  here 
for  brevity. 

We  would  like  to  make  a  few  remarks  before  leaving  this 
section. 

Remark  1:  The  choice  of  Ns  symbols  in  a  group  needs 
to  satisfy  the  condition:  2n  \frj]  NST  <  jg,  to  ensure  that  the 
maximum  rotating  phase  does  not  exceed  a  decision  region 
of  MPSK,  where  |/^|  is  the  absolute  value  of  the  maximum 
possible  Doppler. 

Remark  2:  The  group-wise  estimation  of  the  rotating 
phase  is  insensitive  to  noise  perturbations  due  to  its  averaging 
process  (30),  which  is  an  implicit  low-pass  filtering  process. 

Remark  3:  The  equivalent  composite  Doppler  for  the  M 
channels  at  p- th  group  can  be  approximated  by  f,Pp  =  2  fc. 

VI.  Experimental  Results 

In  this  section,  we  present  results  of  the  proposed  frequency- 
domain  channel  estimation  and  equalization  using  experimen¬ 
tal  data.  The  data  were  collected  off  the  coast  at  Panama  City, 
Florida,  USA,  in  June  2007.  Eight  hydrophones  were  arranged 
unequally  spaced  over  1.86  meters  on  a  vertical  linear  array. 
Both  the  transducer  and  hydrophone  array  were  deployed  at 
the  bottom  of  the  water  with  20  m  depth.  The  source-receiver 
range  is  5.06  km.  QPSK  signals  with  a  bandwidth  of  4  kHz 
(T  =  0.25 ms)  were  transmitted  on  a  carrier  frequency  of  17 
kHz. 


Fig.  1.  Packet  structure.  Each  packet  contains  M  blocks  of  data,  each  data 
block  contains  one  pilot  symbol  and  N—l  information  symbols. 

The  QPSK  signals  were  partitioned  in  packets.  The  packet 
structure  is  depicted  in  Fig.  1.  Each  packet  starts  with  an 


LFM  signal  followed  by  a  gap,  then  a  training  block  signal 
of  511-bit  m-sequence  followed  by  another  gap,  then  the  data 
package  with  M  blocks  of  data,  in  which  each  block  of  data 
consists  of  N  symbols  followed  by  a  guard  time  Tg  seconds, 
finally,  the  packet  ends  up  with  another  LFM  signal  and  a 
gap.  Each  block  of  data  begins  with  one  pilot  symbol  for 
phase  reference.  Therefore,  each  packet  contains  M(N  —  1) 
symbols  of  information  data.  The  gap  after  the  training  signal 
is  long  enough  so  that  the  channel  impulse  response  can  be 
estimated  from  the  training  signal.  The  guard  time  Tg  after 
each  N  symbols  of  data  is  sufficiently  long  for  avoiding  inter 
block  interference.  The  symbol  synchronization  is  carried  out 
by  the  beginning  and  ending  LFM  signals.  The  gap  after  the 
ending  LFM  is  designed  to  avoid  inter-packet  interference. 

The  initial  channel  estimation  was  achieved  using  the  train¬ 
ing  signal  by  the  proposed  frequency-domain  method  detailed 
in  Section  III.  As  an  example,  the  time-domain  amplitude 
responses  of  the  eight  channels  of  the  first  packet  is  depicted  in 
Fig.  2.  The  most  of  the  channel  energy  is  concentrated  within 
10  ms  which  corresponds  to  a  channel  length  of  L  =  40  in 
terms  of  symbol-interval. 


Channel  Impulse  Responses  (8  Channels) 


Time  (ms)  Time  (ms) 


Fig.  2.  A  snapshot  of  the  eight  channel  impulse  responses  of  the  first  packet. 

In  this  paper,  we  report  results  for  the  block  size  of  N  =  512 
and  N  =  1024.  When  the  block  size  N  =  512,  M  was  chosen 
to  be  86,  that  means  each  packet  contains  2M  (N —1)  =  87892 
information  bits.  The  uncoded  average  bit  error  rates  (BERs) 
based  on  4  packets  are  listed  in  Table  1  for  single  channel  and 
multiple  channel  receivers. 

When  the  block  size  N  =  1024,  M  was  chosen  to  be 
48,  that  means  each  packet  contains  2 M(N  —  1)  =  98208 
information  bits.  The  uncoded  average  BERs  based  on  4 
packets  are  listed  in  Table  2  for  single  channel  and  multiple 
channel  receivers. 

As  can  be  seen  from  both  Tables  1  and  2,  the  proposed 
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algorithm  has  better  BER  performance  for  N  =  512  than  that 
of  N  =  1024  when  the  diversity  channel  number  is  less  than 
or  equal  to  4.  However,  when  the  diversity  channel  number  is 
larger  than  or  equal  to  5,  the  proposed  algorithm  has  similar 
BER  performance  for  both  block  size  N  =  512  and  N  = 
1024. 


Table  1:  Average  BER  of  single  channel  and  multiple 
channel  receivers  for  N  =  512. 


Number  of 

Channels  used 

number  of  information 
bits  per  packet 

bit  error 

rate 

1 

87892 

0.48 

2 

87892 

0.075 

3 

87892 

0.013 

4 

87892 

0.0029 

5 

87892 

0.0011 

6 

87892 

4.54  x  10"4 

7 

87892 

7.95  x  10-5 

8 

87892 

6.81  x  10"5 

Table  2:  Average  BER  of  single  channel  and  multiple 
channel  receivers  for  N  =  1024. 


Number  of 

Channels  used 

number  of  information 
bits  per  packet 

bit  error 

rate 

1 

98208 

0.49 

2 

98208 

0.096 

3 

98208 

0.036 

4 

98208 

0.013 

5 

98208 

0.0012 

6 

98208 

2.34  x  10"4 

7 

98208 

9.83  x  10"b 

8 

98208 

8.14  x  10"5 

VII.  Conclusion 

A  novel  frequency-domain  channel  estimation  and  equaliza¬ 
tion  scheme  has  been  proposed  for  single  carrier  underwater 
acoustic  communication  systems.  The  frequency-domain  chan¬ 
nel  estimator  utilizes  a  small-block  training  signal  for  initial 
channel  estimation  and  converts  the  estimated  transfer  function 
matrices  to  a  desired  Appoint  frequency-domain  channel  trans¬ 
fer  function  matrices.  Then  the  frequency-domain  equalizer 
uses  the  linear  MMSE  approach  to  eliminate  the  ISI  caused  by 
the  fading  channel.  The  resulting  signals  are  converted  to  time- 
domain  and  phase  correction  is  performed  with  the  estimated 
group  average  rotating  phases.  The  channel  transfer  function 
matrices  were  re-estimated  by  using  the  detected  data  and 
received  data  of  the  current  block,  these  re-estimated  transfer 
function  matrices  were  employed  to  equalize  the  next  block 
of  data  in  a  decision-directed  manner.  The  proposed  scheme 
has  been  used  to  process  data  measured  by  field  experiments 
conducted  off  the  coast  at  Panama  City,  Florida,  USA,  in  lune 
2007.  In  the  experiment,  QPSK  signals  with  a  bandwidth  of  4 
kHz  were  transmitted  on  a  carrier  frequency  of  17  kHz.  The 


experimental  results  have  shown  that  the  proposed  frequency- 
domain  equalizer  can  achieve  uncoded  BER  on  the  order 
of  10-5  using  7  ~  8  channel  diversity.  It  can  also  achieve 
uncoded  BER  between  1.2  x  10-3  and  4.5  x  10-4  for  5  ~  6 
channel  diversity  systems  under  a  typical  fading  channel  of 
length  10  ms. 
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